METHOD AND SYSTEM FOR SATELLITE BASED PHASE MEASUREMENTS 
FOR RELATIVE POSITIONING OF FIXED OR SLOW MOVING POINTS 
IN CLOSE PROXIMITY 

CROSS-REFERENCE TO RELATED APPLICATIONS 

[0001] This application claims the benefit of U.S. Provisional Application No. 

60/464,756, filed April 23, 2003 April 23. 2003, the contents of which are incorporated by 
reference herein in their entirety. 

BACKGROUND OF THE INVENTION 

[0002] The invention relates generally to Global Positioning System (GPS) receivers and 

more particularly to a method and an apparatus for computing multiple precise locations using 
differential carrier phases of a GPS satellite signal by synchronizing the clocks between the 
master receiver and the slave receiver. It further describes a technique of connecting a plurality 
of antennas to the slave receiver, which can be switched [[in]] onto measure each antennas 
antenna's relative location to the master antenna for monitoring long-term deformation. 

GPS Background 

[0003] The Global Positioning System (GPS) was established by the United States 

government, and employs a constellation of 24 or more satei lites in well-defined orbits at an 
altitude of approximately 26,500 km. These satellites continually transmit microwave L-band 
radio signals in two frequency bands, centered at 1575.42 MHz and 1227.6 MHz., denoted as LI 
and L2 respectively. These signals include timing patterns relative to the satellite's onboard 
precision clock (which is kept synchronized by a ground station) as well as a navigation message 
giving the precise orbital positions of the satellites. GPS receivers process the radio signals, 
computing ranges to the GPS satellites and by triangulating i:hese ranges; the GPS receiver 
determines its position and its internal clock error. Different levels of accuracies can be achieved 
depending on the techniques deployed. This invention specifically targets the sub-centimeter 
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accuracies achievable on a remote and possibly mobile GPS receiver by processing carrier phase 
observations both from the remote receiver and from one or more fixed-position reference 
stations. This procedure is often referred to as Real Time Ki nematic or RTK. 

[0004] To gain a better understanding of the accuracy levels achievable by using the GPS 

system, it is n e c e s s ary necessary to understand the two type;: of signals available from the GPS 
satellites. The first type of signal includes both the Coarse Acquisition (C/A), which modulates 
the LI radio signal and precision (P) code, which modulates, both the LI and L2 radio signals. 
These are pseudorandom digital codes that provide a known pattern that can be compared to the 
receiver's version of that pattern. By measuring the time- shift required to align the 
pseudorandom digital codes, the GPS receiver is able to compute an unambiguous pseudo-range 
to the satellite. Both the C/A and P codes have a relatively long "wav e l e ngth," "wavelength", of 
about 300 meters (1 microsecond) and 30 meters (1/10 microsecond), respectively. 
Consequently, use of the C/A code and the P code yield position data only at a relatively coarse 
level of resolution. 

[0005] The second type of signal utilized for position determination is the carrier signals 

signal . The term "carrier", as used herein, refers to the dominant spectral component which 
remains in the radio signal after the spectral content caused by the modulated pseudorandom 
digital codes (C/A and P) is removed. The LI and L2 carriei signals have wavelengths of about 
19 and 24 centimeters, respectively. The GPS receiver is able to "track" these carrier signals, and 
in doing so, make measurements of the carrier phase to a small fraction of a complete 
wavelength, permitting range measurement to an accuracy of less than a centimeter. 

[0006] In stand-alone GPS systems that determine a receiver's position coordinates 

without reference to a nearby reference receiver, the process of position determination is subject 
to errors from a number of sources. These include errors in the satellite's clock reference, the 
location of the orbiting satellite, ionospheric refraction error:? (which delay GPS code signals but 
advance GPS carrier signals), and tropospheric induced delay errors. Prior to May 2 of 2002 2, 
2002 . a large portion of the satellite's clock error, referred to as Selective Availability (SA) was 
purposefully induced by the U.S. Department of Defense to limit GPS accuracy to non- 
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authorized users. SA would often cause positioning errors exceeding 40 meters, but even today, 
with SA off, errors caused by the ionosphere can be tens of meters. The above mentioned error 
sources (satellite clock and satellite position errors, ionosphere refraction, tropospheric delay and 
SA) are common-mode errors for two receivers that are neaiby. That is, the errors caused by 
these sources are nearly the same for each r e c e iver receiver. 

[0007] Another error source, which is present in the carrier phase measurements, is the 

clock diff e r e nc es difference between the two receivers. This clock difference applies to all 
satellite measurements equally, and as such, can be eliminated by what is known as double 
differencing. This is where one of the satellites is used as a reference and the other satellite 
measurements are compared to it. This reduces the number of usable satellite measurements by 
one. As will be explained later, the more measurements available the better the final solution. 

[0008] To overcome the common-mode errors of the stand-alone GPS system, many 

kinematic positioning applications make use of multiple GPS receivers. A reference receiver 
located at a reference site having known coordinates receives the satellite signals simultaneously 
with the receipt of signals by a remote receiver. Depending on the separation distance, the 
common mod e the common-mode errors mentioned above will affect the satellite signals equally 

for the two receivers. By taking the difference between signals received both at the reference site 

and at the remote location, common-mode errors are effectively eliminated. This facilitates an 

accurate determination of the remote receiver's coordinates relative to the reference receiver's 

coordinates. 

[0009] The technique of differencing signals is known in the art as differential GPS 

(DGPS). The combination of DGPS with precise measurements of carrier phase leads to position 
accuracies of less than one centimeter root-mean- squared (centimeter-level positioning). When 
DGPS positioning utilizing carrier phase is done in real-time while the remote receiver is 
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potentially in motion, it is often referred to as Real-Time Kinematic (RTK) positioning. 

[001 0] One of the difficulties in performing RTK positioning using carrier signals is the 

existence of an inherent ambiguity that arises because each cycle of the carrier signal looks 
exactly alike. Therefore, the range measurement based upon carrier phase has an ambiguity 
equivalent to an integral number of carrier signal wavelengths. Various techniques are used to 
resolve the ambiguity, which usually involves some form of double-differencing of the carrier 
measurements. Once ambiguities are solved, however, the receiver continues to apply a constant 
ambiguity correction to a carrier measurement until loss of lock on that carrier signal or partial 
loss of lock that results in a carrier cycle slip. 

[001 1] Regardless of the technique deployed, the problem of solving integer ambiguities, 

in real-time, is always faster and more robust if there are me re measurements upon which to 
discriminate the true integer ambiguities. Robust means that there is less chance of choosing an 
incorrect set of ambiguities. The degree to which the carrier measurements collectively agree to a 
common location of the GPS receiver is used as a discriminator in choosing the correct set of 
ambiguities. The more carrier phase measurements that are available, the more likely it is that the 
best measure of agreement will correspond to the true (relative to the reference GPS) position of 
the remote GPS receiver. One method, which effectively gives more measurements, is to use 
carrier phase measurements on both L I and L2. The problem though is that it is relatively 
difficult to track L2 because it is modulated only by P code ;ind United States Department of 
Defense has limited access to P code modulation by encrypting the P code prior to transmission. 
Some receivers are capable of applying various cross-correlation techniques to track the P code 
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on L2, but these are usually more expensive receivers [[that]] than LI only capable receivers. 

[0012] Other approaches have been employed to gain additional measurements on GPS 

receivers utilizing additional satellites and other types of satellite systems such as the GLONASS 
system, pseudolites, or Low Earth Orbit (LEO) satellite signals in an attempt to enhance RTK. 
Nevertheless, it is often desired to perform RTK on low-cost LI only receivers that do not have 
access to the GLONASS system, pseudolites, or LEO satellite signals. 

SUMMARY OF THE INVENTION 

[0013] Disclosed herein in an exemplary embodiment is a method for measuring relative 

position of fixed or slow-moving points in close proximity comprising: receiving a set of satellite 
signals with a first receiver corresponding to a first position; receiving a related set of satellite 
signals with a second receiver corresponding to a second position; and computing a position of 
the second position based on at least one of code phase and carrier phase differencing techniques. 
At l e ast one of: Either a clock used in the first receiver and a clock used in the second receiver 
are synchronized to eliminate substantial clock variation between the first receiver and the 
second receiver; and or the first receiver and the second receiver share a common clock. 

[001 4] Also disclosed herein in another exemplary embodiment is a system for 

measuring relative position of fixed or slow-moving points in close proximity comprising: a first 
receiver in operable communication with a first antenna configured to receive a first plurality of 
satellite signals at a first position; and a second receiver in operable communication with a 


second antenna configured to receive a second plurality of satellite signals at a second position; 
and at least one of the first receiver and the second receiver computing a position corresponding 
to a position of the second antenna based on at least one of code phase and carrier phase 
differencing techniques. At least on e of; Either a clock used in the first receiver and a clock used 
in the second receiver are synchronized to eliminate clock variation between the first receiver 
and the second receiver, and or the first receiver and the second receiver share a common clock. 

[0015] Further, disclosed herein in yet another exemplary embodiment is a system for 

measuring relative position of fixed or slow-moving points in close proximity comprising: a 
means for receiving a set of satellite signals with a first receiver corresponding to a first position; 
a means for receiving a related set of satellite signals with a second receiver corresponding to a 
second position; and a means for computing a position of the second position based on at least 
one of code phase and carrier phase differencing techniques. At l e ast on e of: Either a clock used 
in the first receiver and a clock used in the second receiver are synchronized to eliminate clock 
variation between the first receiver and the second receiver, and or the first receiver and the 
second receiver share a common clock. 

[00 1 6] Also disclosed herein in yet another exemplary embodiment is a storage medium 

encoded with a machine-readable computer program code, the code including instructions for 
causing a computer to implement the abovementioned method for measuring relative position of 
fixed or slow-moving points in close proximity. 


[0017] 


Further disclosed herein in yet another exemplary embodiment is a computer data 


signal, the computer data signal comprising code configured to cause a processor to implement 
the abovementioned method for measuring relative position of fixed or slow-moving points in 
close proximity. 
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BRIEF DESCRIPTION OF THE DRAWINGS 

[0018] Referring now to the drawings wherein like elements are numbered alike in the 

several FIGURES: 

[0019] FIG. 1 is a block diagram showing the multiple antennas connected via switches 

to the slave receiver and the single master receiver within the same enclosure to permit clock 
synchronization; 

[0020] FIG. 2 is a diagram depicting signals receive:! from multiple satellites at two 

antenna locations. 

[0020.1] FIG. 3 is a diagram depicting a vessel with multiple antennas and an orientation 
device. 

DETAILED DESCRIPTION 

[0021] This invention discloses the use of two receivers, which either share the same 

clock, or have a clock synchronization technique to eliminate the receiver clock errors. Further 
the reference receiver (herein called the master) is connected to a single antenna whereas the 
slave receiver, which is clock synchronized with the master, has a multitude of antennas 
connected which are switched in and out to take a measurement at each antenna location. 

[0022] The GPS rover receiver computes the location vector from a double or single difference 
of the GPS rover and reference carrier phases for a plurality of GPS satellites. As the receivers 
are either co-located or have a link, the raw m e asur e m e nt measurements from the slave antennas 
are sent to the master for computation (of course any receiver or even a separate computer could 
perform this computation). This eliminates the need for a radio link between the master and slave 
receivers as is required in prior art RTK. 
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[0023] According to a more specific aspect of the present invention, in order to solve the integer 
ambiguity problem, the master selects the slave antenna to be measured based on the GPS 
satellite almanac to provide the best geometry (or one of the best) and based on its time slot. The 
master also has the slave antenna's position stored to provide an immediate calculation of the 
carrier cycle ambiguity to each satellite. Position calculation then follows conventional RTK 
GPS practice of using single or double difference equations involving the total phase distance to 
each satellite to solve the relative location of slave the slave antenna with respect to the master 
antenna. As previously described, there is no clock difference between the two receivers (or the 
clock difference is known and nearly constant) so double differencing may not be required. 
There may however be a significant delay through the coaxial cable to each slave antenna. This 
also can be stored and the delay removed to the measurements. A temperature drift may be 
noticed which will gradually change the delay, but this too can be eliminated by the addition of a 
thermocouple to determine the ambient temperature around the cable and antennas. By doing 
this, all satellite measurements may be used in the solution. 

[0024] Another advantage of eliminating double differencing is that ambiguity search 

routines will not have to form linear combinations to decorrelate the measurement data. When it 
is possible to use single differences, they are generally preferred over double differences 
equations. The double difference cross-correlations are more difficult to deal with 
mathematically, say in a measurement covariance matrix of a Kalman filter. Single difference 
equations result in a measurement covariance matrix having zero cross correlation cross- 
correlation. (But note that if the mathematics [[is]] are handled correctly the accuracy of both 
approaches is the same, it is just that the single difference is easier to handle corr e ctly). 
correctly.) 

[0025] Referring now to FIGS. 1 and 2, a simplified block diagram of the system 10 is 

depicted. In an exemplary embodiment, a method and system to u se of 2 use two receivers, 
which either share the same clock, or include a clock synchronization technique to eliminate the 
receiver clock e rrors is disclosed, errors. Further the reference receiver (hereinafter also called 
the master) 12 is connected to a master antenna, whereas the rover or slave receiver 14, which is 
clock synchronized with the master, has a multitude of antennas 18 connected which are 
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switched in and out to take a measurement at each antenna location. In addition, the mat e r master 
receiver 12 and slave receiver 14 may include direct connection for wireless communication to 
facilitate communication between them. It will be appreciated that while an exemplary 
embodiment is described and illustrated with respect to measuring movement of a dam, dike or 
b e am. The beam, the disclosed invention is readily applicable to other applications where fixed 
or slow moving slow-moving phenomena are tracked. Such applications may include roadways,, 
bridges, building motion, glacier and iceberg travels and the like. It is also applicable to 
conventional RTK applications that require relatively short distance between master and slave 
and where it is desirable to take advantage of a common clock for added robustness and the 
elimination of a radio for cost and robustenss. robustness. For example, one application is local 
surveying or measuring distance at a construction site, or leveling (such as required for 
foundation placement) at that site. 

[0026] In an exemplary embodiment a master receiv er 1 2 also referred to as a reference 

receiver, and a slave receiver 14, also referred to as a rover or remote receiver are substantially 
co-located. The master and slave receivers 12 and 14 respectively, are configured to either share 
the same clock, or include a clock synchronization system. This technique facilitates elimination 
of the receiver clock errors. In an exemplary embodiment, the GPS slave receiver 14 computes a 
location vector based on a double or single difference of the GPS code and/or carrier phases for 
both the master receiver 12 and slave receiver 14 and for a plurality of GPS satellites. As the 
master and slave receivers [[12,]] 12 and 14 are either co-located or have a link, the raw 
measurements from the slave antennas are sent to the master for computation (of course any 
receiver or even a separate computer could perform this computation). This eliminates the need 
for a radio link between the master and slave receivers 12, 14 as is required in existing RTK 
applications. Moreover, in another exemplary embodiment, satellite signals from multiple 
antennas with a known dimensional separation may be combined to achieve receiving an optimal 
set of satellite signals for a given location. Such an approach will be beneficial for instances 
when insufficient data is available from a single antenna or a less desirable set of satellite signals 
[[are]] is all that is available. In this way, a location may still be computed despite poor satellite 
geomet e r , geometry obstructions, and the like. 
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[0027] Advantageously, in an exemplary embodiment, rather than increasing the number 

of measurements, a reduction in the number of unknowns is achieved by eliminating the clock 
errors between the reference receiver 12 and the rover 14 (or master and slave). This approach 
yields an even greater advantage than adding measurements, unless a substantial number of 
measurements could readily be added. In addition, an exemplary embodiment as disclosed herein 
significantly improves the ability to calculate the integer ambiguities to each satellite. [[In]] It 
will be appreciated that because the slave antennas 1 8 are presumed to move far less than a 
fraction of a carrier cycle (e.g., 19 cm) between measuremen ts, the positions of each slave 
antenna 18 location may be stored and then later retrieved an needed to facilitate the immediate 
calculation of the integer ambiguities. 

[0028] In order to solve the integer ambiguity problem with current RTK applications, 

the master receiver 12 selects a particular slave antenna 1 8 to be measured based on the GPS 
satellite almanac to provide the best geometry (or one of the best) and based on its time slot. The 
master receiver 12 also has the slave antenna's position stored (as stated above) to provide an 
immediate calculation of the carrier cycle ambiguity to each satellite. Position calculation then 
follows RTK GPS practice of using single or double difference equations involving the total 
phase distance to each satellite to solve the relative location of slave antenna 1 8 with respect to 
the master antenna 16. One such methodology for GPS positioning employing RTK is taught by 
Whitehead, U.S. Pat. No. 6,469,663 the contents of which are incorporated by reference herein in 
their entirety. As previously described, there is no clock difference between the two receivers 12 
and 14 (or the clock difference is known and nearly constant ) constant), so double differencing 
may not be required. It will however, be readily appreciated chat there may be a significant delay 
through the coaxial cable 20 to each slave antenna 1 8. This delay is dependent upon the selected 
position for each antenna relative to the master (e.g., the length of cable to reach each antenna). 
Advantageously, the delay may readily be measured and stored and the delay mathematically 
removed to correct the measurements. Moreover, selected antennas may exhibit a temperature 
drift [[the]] that may result in a gradual change of the expected delay. However, advantageously, 
this too may be readily eliminated by the addition of a temperature sensor 22 e.g., thermocouple 
and the like, to determine the ambient temperature around the cable 20 and antennas e.g., 16 and 
18. Advantageously, by employing the abovementioned correction and compensation schemes, 
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all satellite measurements may be used to formulate the solution. 

[0029] Another advantage of eliminating double differencing is that ambiguity search 

routines will not have to form linear combinations to decorrelate the measurement data. When it 
is possible to use single differences, they are generally preferred over double differences 
equations. The double difference cross-correlations are more difficult to deal with 
mathematically, say in a measurement covariance matrix of a Kalman filter. Single difference 
equations result in a measurement covariance matrix with zero cross corr e lation cross- 
correlation, which facilitates computation of the ambiguities. It should of course be noted, that if 
the mathematics [[is]] are handled correctly, the accuracy of both approaches is the same. 
However, utilizing the single difference is an easier process. 

[0030] In yet another exemplary embodiment as an enhan cement to the abovementioned 

embodiments, is the capability to take advantage of the slow dynamics of antenna motion by 
averaging over periods of time thereby reducing multipath contributions (which are time 
varying) and poor satellite geometries. In fact, it will be appreciated that the master receiver 12 is 
constantly tracking the satellites and may further be employ id to select the best time of [[day]] 
day, e.g., the best constellation (the GPS satellites orbit in a 12 hour cycle) to perform the 
measurements based on its knowledge of the slave antennas 1 8 position and the satellites 
currently visible. Additionally the master receiver 12 may select two separate times of day[[,]] to 
provide two independent satellite positions position constellations for performing the 
measurements. This would reduce the amount of averaging 1 [me required, yet still provide the 
multipath and poor satellite geometry reduction benefits. Overall, such an approach may be 
employed to reduce power consumption requirements as the receiver would not have to be 
averaging continuously for a twelve hour period. Power consumption reduction is always 
beneficial, especially at remote sites. 

[003 1 ] Referring once again to FIG. 1 , an exemplary embodiment is shown using a 

plurality of slave antennas 18 (also denoted as Al, A2 . . . An) connected to the slave receiver 
14. Each slave antenna 18 is switched (except the last one [[in]] which is selected when all 
switches are connected through to it is selected ) with a switch box 24 (also denoted as SI, S2 . . 
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.). The switch(es) 24 are selected by a controller (in an exemplary embodiment, part of the 
master receiver 12), which may send a tone or some other control signal 30 on the cable 20 to 
activate a particular desired switch 24 and thereby the slave antenna 1 8 connected th e re to 
thereto . It will be appreciated that in order to provide fault protection, the switch(es) 24 may be 
designed and configured so that in the event a switch 24 fails, the connection through to the next 
switch 24 is made. Advantageously, in this way, if one switch 24 should fail, it will still permit 
measurements on the remaining slave antennas 18. As is she wn in the figure, in one exemplary 
embodiment, both the master and the slave receivers 12 and 14 respectively, are integrated on a 
single printed circuit board (PCB), permitting the master and slave receivers to share a common 
clock. Moreover, in an exemplary embodiment, smart reset circuitry is employed to ensure that 
they (the master receiver 12 and slave receiver 14) will start up at exactly the same time and 
therefore the samples will be aligned as well. This approach substantially eliminates the receiver 
clock biases. 

[0032] As mentioned previously, phase drift and delay can result from the coaxial cables, 

which may be removed and/or compensated by using a temperature sensor 22 e.g., a 
thermocouple to measure the temperature. A look-up table may be employed that has stored 
( alternat e ly alternatively a simple formula may be used to save memory) phase delay difference 
versus ambient temperature. An alternative embodiment could use equivalent coaxial cable 
lengths to all antennas including the master so any temperature or other loss and drift effects 
would be matched and therefore cancelled in the single difference calculation. 

[0033] Normally in order to solve for integer ambiguities from [[and]] GPS satellite 

signals, double differencing is used to bring forth the integer nature of the ambiguities by 
removing other non-integer sources of error such as clock and atmospheric delays from the 
measurements. To illustrate, consider four equations describing pseudo-ranges resulting from 
measurements of carrier phase on receivers denoted m and n for the slave and master, 
respectively: 
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K = R m + TSV ' + A + B m + NJ 

<f>j = Rj + tsv' + A 1 + B u + Nj 
(/> m k = Rk + tsv" = Ak + B m + N m k 

(j>* = R n k + rsv k = A k + Bn + N n k [[i.]] 

[0034] Here $ m ' is the measured pseudorange from rover receiver m to satellite i, <f)^ is 

the measured pseudorange from reference receiver n to satel lite /', (j) m is the measured 
pseudorange from rover receiver m to satellite k, and $ n is tlie measured pseudorange from 
reference receiver n to satellite k. Each pseudorange is actually a measure of the summation of a 
number of different physical quantities all of which shall be expressed in units of carrier cycles at 
LI (roughly 19 cm). 

[0035] Specifically, in the first of these equations, the term Rj is the true geometric 

range from receiver m to satellite i, tsv' is the clock error of satellite /, A' is the atmospheric 
delays, which are associated with satellite z, B m is the clock error of receiver m, and Nj is the 
integer ambiguity in the range measurement from receiver m to satellite /. Similar notation 
applies to the remaining three equations. For simplicity, these equations do not show noise 
effects such as errors caused by receiver thermal noise or multipath noise. 

[0036] Consider first applying the single difference. 1 f the first two equations are 

differenced: 

/. ^ - (/) ' = R j - Rj + B m - B n + N,j - Nj 

ii. Similarly, differencing the second two equations yields: 
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iii. <j) k - <j>* = R m k - R k +B m - B n + N m k - N n k 

[0037] The satellite common errors, such as satellite clock, isv' and atmosphere, A' 

(atmosphere is common if we assume relative close proximity of receivers m and n) are removed 
in the single difference. As the clock errors B m are common these term will also cancel out, 
leaving: 

^ - = Rj - r 1 + nJ 

[0038] Since the ambiguities are all integers that can be lumped together into a single 

term, it may be written: 

- - Rj - Rj + N mn 

where 

N mn = NjNj 

[0039] This shows that single differencing the pseudorange measurements removes 

common atmospheric errors from the equations while leaving simple combinations of the 
geometric ranges and integer ambiguities, and clock errors-d-eps drop out due to the 
synchronization of the two receivers. For N satellites in common view of the master (reference) 
and slave (remote) receivers 12 and 14 respectively, there are N such single-difference equations 
that can be formed without causing mathematical redundancy. Whereas double differencing, to 
eliminate clock biases in receivers, which are not clock synchronous, results in only N-l 
equations. This gives rise to N unknown integer ambiguities that must be solved in addition to 
the [[3]] three u nknown coordinates (X, Y, Z) of the GPS receiver. Note that each geometric 
range term, for example R' m , is a function only of the receiver's position and the transmitting 
satellite's position. Specifically: 
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^{Xrecv m - Xsat') 2 + (Yrecv m - Ysat'f + (Zrecv m - Zsat 1 ) 2 


[0040] where Xrecv m , Yrecv m Zrecv m are the Cartesian coordinates of the receiver m at 

the time reception of the signal from satellite i, whose coordinates are Xsat' , Ysat' , Zsat' at 
the time of signal transmission. In the problem at hand, only the selected slave's antenna's 18 
position is unknown. Once the ambiguities are determined, only the selected antenna's 3- 
coordinates of position are unknown and these are easily solved using a mathematical approach 
such as Least Squares. 

[0041] Every time a new slave antenna 1 8 is selected, the integer ambiguities must be 

solved. This is a complex process and can be very time consuming if the position is unknown. 
However, in this instance, it will be appreciated that the movements to be measured are on the 
order of less than a quarter of a wavelength (5 cm) between measurements. This limitation 
permits a rapid calculation of the integer ambiguities since the master receiver 12 "knows" the 
satellite's position and the selected antenna's position well enough to directly calculate 
ambiguities. Such an approach will greatly reduce the time utilized to solve for the integer from 
up to 10 minutes to a second or less. Cycle slips, which result usually from motion which the 
receiver failed to track properly and therefore slipped from ocie ambiguity to another is also 
greatly reduced due to the very low dynamics of the selected antenna location. An added benefit 
of the low dynamics is the receiver can integrate the measurements over a long period of time 
and narrow the carrier tracking loop bandwidth to reduce noise. 

[0042] As mentioned previously, it should be appreciated that another source of error in 

applying RTK positioning, especially when solving for integer ambiguities over long baselines, 
is non-common atmospheric propagation delays on the signals received by the slave (rover) 14 
and master (reference) receivers 12. Since differencing cannot eliminate these non-common 
delays, the next best alternative is to estimate or model their effects. However, [[In]] in an 
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exemplary embodiment, the slave antennas 18 and the master antenna 16 will, most likely, be 
within 5 kilometers of each other and at this distance the atmospheric effects are minimal and 
may readily be ignored. 

[0043] A further advantage of this technique should permit a carrier phase based phase- 

based solution even when a large portion of the sky, and therefore the visible satellites, are 
obscured by a wall, dam or other structure. This is because, as described above, the receiver will 
still have one more measurement than previously due to the utilization of single differencing 
rather than double differencing t e chnique techniques . In addition, the fixed or very slow moving 
nature of the problem permits long-term measurements. 

[0044] Referring now to FIG. 2 as well, in yet another exemplary embodiment, a 

technique is employed to utilize and take advantage of the master receiver's 12 knowledge of the 
satellite's location in the sky, and a preprogrammed knowledge of the visibility of the sky for 
selected slave antennas 1 8. The master receiver 12 may then [chose] choose the best time, that is, 
the time with the most satellites visible to the selected slave antenna 18, to perform the 
measurement at that location. The receiver can then dwell for some time (say one half hour) to 
integrate and reduce noise, then move on to another slave antenna 1 8. Moreover, it will be 
appreciated that the master receiver 12 may direct that the slave receiver return to the same 
location after some duration ergr e.g., a few hours, when another optimal/desirable geometry is 
available, which is uncorrelated to the first. By taking measurements at two (or more) different 
times (and geometries), and averaging the two (or more) measurements, multipath and 
atmospheric induced errors, typically correlated over time, will be reduced. This method will 
allow monitoring of the face of a dam or berm, or even a valley wall, which was previously 
impossible to monitor. 

[0045] Further assumptions may be made of the anticipated motion of the monitoring 

point at the selected slave antenna 1 8 to further reduce the number of measurements required. 

[0046] For example, if it is a dam, the anticipated moi:ion is horizontally away from the 

pressure excerpted by the material behind the dam. By performing the calculation only [[on]] in 
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this direction, a single satellite may be enough to perform a measurement. This is obvious when 
looking at this equation: 



y[{Xrecv m - Xsat'f + {Yrecv m - Ysat') 2 + (Zrecv m - Zsat ] ) 2 

[0047] As explained previously the satellite position (Xsat, Ysat and Zsat) are known, 

and if the receiver assumes there is minimal motion in Y arcing and Z, then there is only one 
unknown left. Of course, additional satellites are highly desired to reduce noise and errors and to 
help detect any false or erroneous readings from throwing the solution off. 

[0048] Another area of concern for running a long length of coaxial cable 20 to the 

antennas 16, 18, other than phase delay, which was addressee! earlier, is attenuation. In yet 
another exemplary embodiment, the slave antennas 1 8 may be configured as active antennas, 
e.g., antennas that include an internal Low Noise Amplifier (LNA). In a receiver design, the 
Noise Figure noise figure is often important, and comprises £. combination of the noise 
temperature before the first LNA, the LNA Nois e Figur e noi: ;e figure and subsequent losses 
divided by the LNA gain. Subsequent amplifier gains will reduce following noise temperature 
(T) contributions by their gain as is shown in the equation below: 

Tt=T(preLNA) + T(LNA)+T(lna2)/(CL x.Glnal) +T(lna3)/(CL x Glnal x 
Glna2)+T(lna4)/(CL x Glnal x Glna2 x Glna3) etc. 

[0049] where: CL refers to cable losses in linear terms, that is -10 dB is 0.1, 

[0050] Glnan refers to gain of LNAn in linear terms :;o a gain of 20 dB is 100, 

[005 1 ] T(LNAn) refers to the noise temperature in Kelvin of stage n. 


[0052] 


Noise Figure (F) is related to noise temperature by: 


49- 


F(dB)=10 x L0G((1+T)/Tamb) 

[0053] Where Tamb refers to the reference temperatu re, typically 290 K (20 Celsius). 

[0054] As an example, a typical low loss coaxial cable (RG6 type) has 20 dB (CL=0.01) 

of attenuation every 100 meters. The noise temperature of the antenna and LNA is 170 K (2 dB 
noise figure), the gain of the first LNA is 30 dB (or 1000). Subsequent LNA's have the same 
noise temperature and a gain of 12 dB (15.8). If each antenna is 50 meters apart the losses are - 
10 dB. After five stages the noise temperature of the system is: 

T5=T1+T2/(CL1 x.Gl) + T3/(CL1 x C12 x Gl x.G2)+ T4/(CL1 x CL2 x C13 x 
Gl xG2 G3) + T5/(CL1 xC12xC13xC14xGl xG2xG3xG4) 

T5= 1 90+ 1 90/ 1 00+ 1 90/ 1 5 8+ 1 90/250+ 1 90/395 

T5=194K 

F5=2.22 dB 

[0055] This is compared to the first stage, which woul d have a noise figure of 2 dB. A 

GPS receiver such as the master receiver 12, or slave receiver 14 can operate with a noise figure 
of up to 3.5 dB without suffering significant degradation. As can be seen, additional stages will 
have diminishing contributions. The total gain will be increasing by only 2 dB each step, so after 
1 km, in this example, the maximum gain will be 68 dB, the gain of the first stage is 30 dB, the 
Automatic Gain Control of the receiver can remove this difference easily. Also after 20 stages (1 
km) the total noise temperature in this example would be T(l km)=T94.7 K, an insignificant 
increase. 


[0056] Further, in another exemplary embodiment, multiple antennas can be used to 

compute a solution of a single point on a rigid body to which they are attached, using known 
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geometry and distances. Such an approach may be employed, for example, when not any one 
antenna provides enough useful information (satellites) to compute a location solution due to 
obstructions, but the conglomerate could. Advantageously, a position solution employing this 
approach would not necessarily have to utilize carri e r phas e b ased carrier phase-based differencing 
(it could be code phase). An application might include positioning on a barge 30 shown in Fig. 3, 
where location is needed but there are many cranes and towers blocking the view so that there is not 
one optimum GPS location. However, by placing antennas Al and A2 respectively on either side of 
the barge 30, enough satellites could be tracked by the combined antenna arrangement that a 
solution of the location of some point on the barge could still be obtained. Furthermore, on the 
barge 30, a compass 32 can also be used to give orientation, thus removing another unknown from 
the relative location of the two antennas. Rather than solving a relative location of one antenna with 
respect to another, the combined outputs of either antenna Al or A2 and the compass 32 can be 
used to produce one non-relative location. 

[0057] It will be appreciated that the satellite systems as discussed herein may include 

but not be limited to Wide Area Augmentation System (WAAS), Global Navigation Satellite 
System (GNSS) including GPS, GLONASS and other satellite ranging technologies. The term 
WAAS here is used herein as a generic reference to all GNSS augmentation systems which, to 
date, include three programs: WAAS (Wide Area Augmentation System) in the USA, EGNOS 
(European Geostationary Navigation Overlay System) in Europe and MSAS (Multifunctional 
Transport Satellite Space-based Augmentation System) in Japan. Each of these three systems, 
which are all compatible, consists of a ground network for observing the GPS constellation, and 
one or more geostationary satellites. 

[0058] It will be appreciate s appreciated that while a particular series of steps or 

procedures is described as part of the abovementioned process, no order of steps should 
necessarily be inferred from the order of presentation. For example, the process includes 
receiving one or more sets of satellite signals. It should be evident the order of receiving the 
satellite signals is variable and could be reversed without impacting the methodology disclosed 
herein or the scope of the claims. 
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[0059] It should further be appreciated that while an exemplary partitioning functionality 

has been provid e d, provided, [[It]] it should be apparent to one skilled in the art [[,]] that the 
partitioning could be different. For example, the control of the master receiver 12 and slave 
receiver 14, could be integrated in any, or another unit. The processes may, for ease of 
implementation, be integrated into a single unit. Such configuration variances should be 
considered equivalent and within the scope of the disclosure and claims herein. 

[0060] The disclosed invention may be embodied in the form of computer-implemented 

processes and apparatuses for practicing those processes. The present invention can also be 
embodied in the form of computer program code containing instructions embodied in tangible 
media, such as floppy diskettes, CD-ROMs, hard drives, or any other computer-readable storage 
medium, wherein, when the computer program code is loaded into and executed by a computer, 
the computer becomes an apparatus for practicing the invention. The present invention can also 
be embodied in the form of computer program code, for example, whether stored in a storage 
medium, loaded into and/or executed by a computer, or as data signal transmitted whether a 
modulated carrier wave or not, over some transmission medium, such as over electrical wiring or 
cabling, through fiber optics, or via electromagnetic radiation, wherein, when the computer 
program code is loaded into and executed by a computer, the computer becomes an apparatus for 
practicing the invention. When implemented on a general-puipose microprocessor, the computer 
program code segments configure the microprocessor to create specific logic circuits. 

[0061] While the description has been made with reference to exemplary embodiments, it 

will be understood by those of ordinary skill in the pertinent art that various changes may be 
made and equivalents may be substituted for the elements thereof without departing from the 
scope of the disclosure. In addition, numerous modifications may be made to adapt the teachings 
of the disclosure to a particular object or situation without departing from the essential scope 
thereof. Therefore, it is intended that the Claims not be limited to the particular embodiments 
disclosed as the currently preferred best modes contemplated for carrying out the teachings 
herein, but that the Claims shall cover all embodiments falling within the true scope and spirit of 
the disclosure. 


